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[i]  The  atmosphere’s  response  to  the  total  solar  eclipse  of  4  December  2002  is  studied 
using  a  prototype  high-altitude  global  numerical  weather  prediction  model 
(NOGAPS-ALPHA).  Local  reductions  in  solar  ultraviolet  (UV)  radiation  during  the 
eclipse  are  estimated  using  astronomical  calculations  of  umbral  and  penumbral  surface 
trajectories  and  observed  solar  limb  darkening  at  ~200-300  nm.  In  NOGAPS-ALPHA 
these  UV  eclipse  shadows  yield  stratospheric  radiative  cooling  rate  footprints  peaking  near 
27  K  day-1,  a  value  2-3  times  larger  than  assumed  in  previous  modeling.  Difference 
fields  between  NOGAPS-ALPHA  runs  with  and  without  this  eclipse  forcing  reveal 
vertically  deep  middle  atmospheric  responses,  with  three-dimensional  horizontal 
structures  very  similar  to  the  large-scale  “bow-wave”  response  first  proposed  by 
Chimonas  (1970).  Such  structure  appears  clearly  only  at  later  times  when  total  eclipses 
have  abated  and  gravity  waves  generated  in  the  stratosphere  have  had  time  to  propagate 
vertically.  Bow-wave  amplitudes  and  direct  thermal  cooling  responses  are  both  small 
(;S1  K  for  temperature  and  ;S2-3  m  s  1  for  horizontal  winds),  contradicting  some 
rocketsonde  measurements  that  suggest  much  larger  responses  near  50-60  km  altitude. 

We  also  find  clear  evidence  of  a  bow-wave-like  response  in  the  model’s  surface 
pressure  fields,  with  an  amplitude  ~0. 1-0.5  hPa,  while  surface  air  temperatures  in 
NOGAPS-ALPHA  show  ~4  K  cooling  over  Africa  during  the  eclipse.  Both  findings  are 
consistent  with  surface  atmospheric  data  acquired  during  previous  eclipse  passages. 
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the  total  solar  eclipse  of  4  December  2002  simulated  with  a  high-altitude  global  model,  J ■  Geophys.  Res.,  112,  D14105, 
doi:  1 0. 1 029/2006 JD007880. 


1.  Introduction 

[2]  Solar  eclipses  provide  natural  predictable  perturbation 
experiments  for  studying  the  atmosphere’s  response  to 
changes  in  solar  forcing.  From  the  earliest  days  of  atmo¬ 
spheric  research,  meteorological  measurements  have  been 
stationed  within  predicted  paths  of  eclipse  shadows  to 
investigate  these  sensitivities  (for  a  brief  historical  review, 
see  A  pi  in  and  Harrison  [2003,  section  2]).  In  the  middle 
atmosphere,  the  earliest  such  measurements  focused  on 
ozone  changes  during  eclipses,  given  evolving  understand¬ 
ing  of  the  role  of  solar  ultraviolet  photolysis  in  odd  oxygen 
production  in  the  stratosphere  and  mesosphere  [Kawabata, 
1937;  Jerlov  et  al.,  1954;  Stranz,  1961;  Hunt,  1965; 
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Randhawa,  1968].  Rocketsonde  sounding  experiments  also 
searched  for  eclipse-related  changes  to  stratospheric  and 
mesospheric  winds  and  temperatures  [ Ballard  et  al,  1969]. 

[3]  These  exploratory  measurements  were  given  greater 
impetus  and  focus  by  the  theoretical  study  of  Chimonas 
[1970],  Modeling  the  eclipse  shadow  as  a  moving  three- 
dimensional  cooling  rate  perturbation  to  the  stratosphere 
due  to  reduced  shortwave  ozone  heating,  Chimonas  [1970] 
solved  the  perturbation  fluid  equations  to  infer  the  atmos¬ 
phere’s  response  to  this  diabatic  forcing.  His  solutions 
predicted  a  long-wavelength  three-dimensional  bow-wave¬ 
like  gravity  wave  in  the  wake  of  the  moving  eclipse  shadow, 
which  propagated  vertically  to  both  higher  and  lower 
altitudes. 

[4]  Just  prior  to  the  total  solar  eclipse  of  7  March  1970, 
Chimonas  and  Hines  [1970]  used  the  predicted  path  of  the 
surface  eclipse  shadow  and  fonnulas  from  Chimonas  [1970] 
to  predict  geographical  regions  where  the  hypothetical  bow 
wave  response  might  be  observed,  both  near  the  ground  and 
in  the  ionosphere.  After  the  event,  Chimonas  and  Hines 
[1971]  argued  that  traveling  ionospheric  disturbances 
(TIDs)  measured  remotely  by  Davis  and  da  Rosa  [1970] 
from  a  ground  location  in  the  path  of  the  eclipse  shadow 
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directly  validated  their  predictions.  However,  more  detailed 
postanalyses  of  the  TIDs  observed  by  Davis  and  da  Rosa 
[1970]  and  others  led  to  a  consensus  that  the  source  could 
not  be  definitively  ascribed  to  the  eclipse:  one  of  the  many 
other  sources  of  TIDs,  such  as  the  strong  magnetic  storm 
activity  that  occurred  at  the  time,  could  just  as  easily  have 
generated  these  waves  [see,  e.g.,  Sears,  1972;  Arendt,  1972; 
Schodel  et  al.,  1973;  Frost  and  Clark,  1973],  Near  the 
surface,  Anderson  et  al.  [1972]  measured  atmospheric 
pressure  oscillations  qualitatively  similar  to  those  predicted 
by  Chimonas  and  Hines  [1970],  but  with  amplitudes  1-2 
orders  of  magnitude  larger.  They  reported  similar  features  in 
other  surface  pressure  measurements  during  eclipses, 
extending  back  to  the  turn  of  the  century.  This  led  Chimonas 
[1973]  to  propose  a  completely  different  theory  for  these 
observations,  in  terms  of  a  Lamb  wave  driven  by  net 
radiative  cooling  of  tropospheric  cloud  layers.  Recent  mod¬ 
eling  has  suggested  that  in  situ  generation  of  atmospheric 
waves  by  eclipses  can  also  occur  in  the  thermosphere 
[Miiller-Wodarg  et  al,  1998],  because  of  reduced  heating 
by  absorption  of  extreme  ultraviolet  solar  radiation. 

[5]  These  source  ambiguities  have  complicated  subse¬ 
quent  experimental  studies  that  have  sought  evidence  of 
eclipse-generated  gravity  waves.  In  the  upper  atmosphere, 
some  observers  report  null  detections  [e.g.,  Schodel  et  al., 
1973;  Hunter  et  al.,  1974;  Boitman  et  al.,  1999],  whereas 
others  find  oscillations  apparently  consistent  with  an 
eclipse-generated  gravity  wave,  but  disagree  on  whether 
the  source  is  stratospheric  or  thermospheric  [e.g.,  Bertin  et 
al.,  1977;  Hanuise  et  al. ,  1982;  Altadil  et  al.,  2001;  Farges 
et  al,  2003;  Jones  et  al.,  2004].  Similarly,  some  studies  of 
near-surface  pressure  observations  present  evidence  of 
eclipse-related  wave  oscillations,  but  disagree  as  to  the 
eclipse-mediated  source  mechanism  [e.g.,  Goodwin  and 
Hobson,  1978;  McIntosh  and  Revelle,  1984;  Seykora  et 
al,  1985;  Farges  et  al.,  2003],  while  others  report  null 
detections  [e.g.,  Jones  and  Bogart,  1975;  Anderson  and 
Keefer,  1975;  Jones,  1976,  1999].  On  collating  results  from 
those  studies  reporting  positive  wave  detections,  both  in  the 
lower  and  upper  atmospheres,  a  wide  spread  of  inferred 
wave  properties  is  revealed  [see,  e.g.,  Farges  et  al.,  2003]. 
These  facts  led  Davies  [1982]  to  question  whether  any 
definitive  experimental  evidence  exists  for  a  characteristic 
eclipse-generated  gravity  wave  due  to  reduced  ozone  heat¬ 
ing  in  the  stratosphere. 

[6]  Direct  upper  stratospheric  or  mesospheric  measure¬ 
ments  could  help  to  clarify  this  situation,  but  far  fewer 
measurements  exist  at  these  altitudes  because  of  the  diffi¬ 
culties  in  acquiring  high-resolution  day-night  data.  Some 
limited  rocketsonde  measurements  have  shown  surprisingly 
large  temperature  decreases  and  associated  meridional  wind 
changes  in  the  lower  mesosphere  during  eclipses  [ Ballard  et 
al.,  1969;  Quiroz  and  Henry,  1973;  Randhawa,  1974; 
Schmidlin  and  Olsen,  1984],  whereas  other  measurements 
have  found  no  discemable  changes  [Randhawa,  1973;  Ball 
et  al.,  1980]. 

[7]  In  summary,  despite  decades  of  research,  observational 
evidence  for  a  characteristic  bow-wave  response  of  the 
atmosphere  to  eclipse  passages  remains  equivocal.  One 
issue  could  be  the  simplicity  of  existing  models.  For 
example,  the  only  extension  to  the  analytical  Chimonas 
[1970]  theory  was  provided  by  Fritts  and  Luo  [1993],  who 


derived  numerical  solutions  to  the  pseudo-incompressible 
perturbation  equations  using  approximations  for  thermal 
eclipse  forcing,  shadow  speed  and  background  atmospheric 
conditions  similar  to  those  used  by  Chimonas  [1970],  Jones 
[1999]  has  argued  for  more  realistic  model  predictions  for 
specific  eclipse  events  to  help  guide  future  observations. 
Several  thermospheric  general  circulation  models  (GCMs) 
have  simulated  the  upper  atmosphere’s  response  to  eclipses, 
finding  strong  in  situ  temperature,  wind  and  composition 
changes  peaking  at  ~250-300  km  altitude  [Ridley  et  al., 
1984;  Roble  et  al.,  1986;  Miiller-Wodarg  et  al.,  1998].  Since 
their  lower  boundaries  were  all  >80  km,  these  models  did 
not  simulate  stratospheric  ozone  cooling  due  to  the  eclipse 
and  thus  could  not  simulate  the  specific  bow  waves  pre¬ 
dicted  by  Chimonas  [1970]  and  Fritts  and  Luo  [1993]. 
Several  mesoscale  weather  models  have  simulated  regional 
changes  in  near-surface  conditions  associated  with  eclipse 
passages  [Gross  and  Hense,  1999;  Prenosil,  2000;  Vogel  et 
al.,  2001],  These  models  do  not  extend  through  the  strato¬ 
sphere  and  lack  the  global  scale  needed  to  simulate  the  fast- 
moving  large-scale  gravity  wave  and  Lamb  wave  responses 
anticipated  by  Chimonas  [1970,  1973], 

[8]  In  this  paper,  we  use  a  state-of-the-art  global  numer¬ 
ical  weather  prediction  (NWP)  model  that,  as  described  in 
section  3,  contains  the  necessary  physics  to  simulate  real¬ 
istic  time-dependent  thermal  cooling  of  the  troposphere  and 
stratosphere  during  a  solar  eclipse.  In  section  2  we  use 
precomputed  time  histories  of  the  lunar  shadow’s  motion 
across  the  Southern  Hemisphere  during  the  total  solar 
eclipse  of  4  December  2002  to  specify  corresponding 
reductions  in  solar  ultraviolet  insolation  during  the  eclipse. 
In  section  4,  difference  fields  between  model  hindcasts  run 
with  and  without  these  eclipse  effects  included  are  used  to 
isolate  the  model’s  simulated  response  to  this  eclipse.  In 
section  5  we  compare  these  responses  to  the  predictions  of 
stratospheric  bow  wave  theory,  tropospheric  Lamb  wave 
theory,  and  relevant  observations.  Major  findings  are 
summarized  in  section  6. 

2.  Total  Solar  Eclipse  of  4  December  2002 

[9]  Detailed  calculations  of  solar  eclipse  properties  are 
issued  ahead  of  time  by  the  U.S.  Naval  Observatory  (USNO). 
Figure  la  characterizes  properties  of  the  4  December  2002 
total  eclipse  from  USNO  calculations.  Gray  curves  show 
30  s  time  series  of  the  southern  and  northern  limits  of  the 
penumbral  (partial)  and  umbral  (total)  eclipse  shadows,  the 
latter  showing  the  very  small  geographical  size  of  the  total 
eclipse  regions  at  any  given  time. 

[10]  Standard  USNO  eclipse  calculations  also  include 
30  min  time  series  of  the  leading  and  lagging  limbs  of  the 
lunar  shadow  on  the  Earth,  which  define  the  geographical 
perimeter  of  the  eclipse  penumbra.  For  this  study,  those 
calculations  were  repeated  to  increase  the  density  of  those 
fields  to  1  min.  However,  these  calculations  do  not  currently 
provide  the  geographical  distribution  of  eclipse  magnitudes 
Em  within  the  penumbral  region.  Appendix  A  describes  the 
calculations  we  performed  to  estimate  time-varying  eclipse 
magnitudes  EM  over  the  Earth’s  surface,  and  to  convert 
them  into  a  corresponding  reduction  factor  for  the  total 
incoming  UV  solar  irradiance,  Euv,  hereafter  called  the  UV 
obscuration. 
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(a)  4  Dec  2002  Solar  Eclipse 


(b)  05:00  UTC 


(f)  07:00  UTC 


(j)  09:00  UTC 


(c)  05:30  UTC 


(g)  07:30  UTC 


(k)  09:30  UTC 


(d)  06:00  UTC 


(h)  08:00  UTC 


(I)  10:00  UTC 


Figure  1.  (a)  Properties  of  the  4  December  2002  solar  eclipse.  Thick  gray  curves  show  northernmost 
and  southernmost  limits  of  eclipse  penumbra  (partial  eclipse),  thin  gray  curve  shows  ground  track  of  the 
eclipse  umbra  (total  eclipse).  Solid  circles  show  points  every  hour  along  these  curves  from  0600  UTC  to 
0900  UTC.  (b-1)  White  contours  show  eclipse  UV  obscurations,  Euv,  at  30  min  intervals  from  0500  UTC 
to  1000  UTC.  Contour  levels  are  0.01,  0.1,  0.3,  0.5,  0.7  and  0.9.  Gray  shading  shows  cosines  of  local 
solar  zenith  angles  (f>  at  cos  cf>  intervals  of  0.1.  The  terminator  (cos  <f>  =  0)  is  labeled  in  each  plot.  A  color 
animation  of  the  time  sequence  in  Figures  lb- 11  is  provided  as  auxiliary  material. 


[ i i ]  Figures  lb- 11  plot  these  USNO-based  Euv 
estimates  every  30  min  during  the  eclipse,  along  with 
contours  of  cos  (j),  where  6  is  the  solar  zenith  angle. 
Animation  S 1 1  provides  a  color  version  of  this  same 
presentation  at  2  min  resolution  for  the  entire  period  of 
the  eclipse,  depicting  both  the  time  evolution  and  speed  of 
passage  of  these  Euv  contours  across  the  globe.  At 
0500  UTC  and  0530  UTC,  partial  eclipsing  of  the  Sun 
commences  over  Africa  during  early  morning  hours,  and  by 
0600  UTC  regions  of  totality  form  over  parts  of  Africa. 
From  0700  to  0800  UTC  the  penumbral  perimeter  closes 
and  the  total  eclipse  regions  move  rapidly  eastward  across 
the  Southern  Ocean  at  local  times  near  midday  (cos  <f>  ~  1). 
Total  eclipses  persist  to  0900  UTC,  then  the  remaining 
partial  eclipse  regions  shrink  and  eventually  disappear  with 
the  setting  Sun  over  Australia  just  after  1000  UTC. 


1  Auxiliary  materials  are  available  in  the  HTML.  doi:10.1029/ 
2006JD007880. 


[12]  Figure  2  plots  the  speed  of  the  umbral  shadow,  V, 
across  the  surface  of  the  Earth.  The  umbral  footprint  moves 
very  rapidly  to  the  east  at  the  start  and  end  times  of  total 
eclipse  in  regions  where  local  solar  zenith  angles  are  large 
and  the  Sun  is  low  in  the  sky  (see,  e.g.,  Figure  Id).  Umbral 
motion  reaches  a  minimum  speed  of  ~70Q  m  s-1  near 
0730  UTC  over  the  Southern  Ocean  between  Africa  and 
Australia  when  the  Sun  is  much  higher  in  the  sky  (see 
Figure  lg). 

3.  Modeling  Tools 
3.1.  NOGAPS-ALPHA 

[13]  The  Navy  Operational  Global  Atmospheric  Prediction 
System  (NOGAPS)  is  the  Department  of  Defense’s  (DoD) 
global  numerical  weather  prediction  system  [ Hogan  and 
Rosmond ,  1991].  FI  ere  we  use  a  prototype  advanced-level 
physics  high-altitude  (ALPFIA)  version  of  the  NOGAPS 
global  spectral  forecast  model  component,  referred  to  as 


3  of  22 


D14105 


ECKERMANN  ET  AL.:  MODELING  ATMOSPHERIC  EFFECTS  OF  A  SOLAR  ECLIPSE 


D14105 


4  December  2002 


0600  0630  0700  0730  0800  0830  0900 
time  of  total  eclipse  (UTC) 


Figure  2.  Speed  V  of  the  umbral  shadow  as  a  function  of 
universal  time  for  the  total  solar  eclipse  of  4  December 
2002. 


NOGAPS-ALPHA.  Complete  descriptions  of  NOGAPS- 
ALPHA  are  provided  by  Eckermann  et  al.  [2004]  and  Allen 
et  al.  [2006],  We  focus  here  only  on  those  components 
salient  to  the  eclipse  simulations  reported  here. 

3.1.1.  Radiation 

3. 1.1.1.  Standard  Model 

[14]  NOGAPS-ALPHA  currently  uses  the  so-called 
CL1RAD  (climate  radiation)  radiative  heating  and  cooling 
rate  parameterizations  of  Chou  and  Suarez  [1999]  and  Chou 
et  al.  [2001],  respectively.  As  implemented  in  NOGAPS- 
ALPHA,  the  Chou  and  Suarez  [1999]  scheme  calculates  the 
absorption  of  ultraviolet,  visible  and  infrared  solar  radiation 
by  03,  02,  C02,  water  vapor  and  both  convective  and  stable 
clouds,  as  well  as  Rayleigh,  surface  and  cloud  scattering.  At 
the  surface,  model  albedos  are  used  to  compute  the  direct 
and  diffuse  infrared,  visible  and  ultraviolet  fluxes.  We  use  a 
solar  constant  So  with  a  base  value  of  1365  W  m~2  that  then 
varies  with  the  Earth’s  orbital  position  according  to  Paltridge 
and  Platt  [1976,  equation  (3.3)].  The  Chou  et  al.  [2001] 
longwave  cooling  scheme  computes  transmission  and  absorp¬ 
tion  by  03,  C02,  water  vapor  and  clouds,  as  well  as  cloud 
scattering.  Both  schemes  yield  heating  and  cooling  rates 
accurate  to  within  5%  of  line-by-line  calculations  from  the 
ground  to  0.01  hPa. 

[15]  The  radiation  schemes  use  NOGAPS-ALPHA’s  prog¬ 
nostic  specific  humidity  and  cloud  fields  between  the  surface 
and  1 00  hPa.  At  altitudes  above  1 00  hPa,  specific  humidities 
are  set  in  the  radiation  schemes  by  zonal-mean  climatologies 
that  vary  as  a  function  of  month,  latitude  and  pressure,  which 
are  interpolated  to  the  day  of  year  and  location  of  the  grid 
point  in  question.  The  climatology  is  based  on  multiyear 
version  18  Halogen  Occultation  Experiment  (HALOE)  data 
[Jackson  et  al.,  1998]  from  100  to  0.3  hPa,  and  diumally 
averaged  output  from  the  CHEM2D  model  [McCormack 
and  Siskind,  2002]  from  0.3  to  0.001  hPa. 

[16]  While  NOGAPS-ALPHA  has  a  prognostic  ozone 
capability  [McCormack  et  al.,  2004],  for  the  runs  reported 
here  we  specified  ozone  for  the  radiation  calculations  using 


a  zonal-mean  observational  climatology  that  varies  with 
month,  latitude  and  pressure,  which  is  then  interpolated  to 
the  specific  day  of  year,  pressure  and  latitude  of  each  model 
grid  box.  The  climatology  used  here  is  an  amalgam  from 
various  sources:  (1)  the  Fortum  and  Kelder  [1998]  clima¬ 
tology  (1000-200  hPa),  (2)  the  1 1-year  HALOE  climatology 
of  Grooss  and  Russell  [2005]  (150-0.1  hPa),  and  (3)  zonal- 
mean  values  from  the  four-dimensional  ozone  climatology 
of  the  United  Kingdom  (UK)  Universities  Global  Atmo¬ 
spheric  Modelling  Project  (UGAMP)  [Li  and  Shine,  1995] 
(0.1-0.001  hPa).  The  resulting  ozone  mixing  ratios  for 
December  are  plotted  in  Figure  3a.  The  values  above  the 
1  hPa  level  are  daytime  averages  only,  mainly  data  from  the 
Solar  Mesosphere  Explorer  in  the  case  of  the  UGAMP 
climatology  [Li  and  Shine,  1995],  and  compare  favorably 
with  more  recent  satellite  measurements  of  daytime  ozone 
in  the  upper  stratosphere  and  mesosphere  [e.g.,  Ricard  et  al., 
1996;  Kaufmann  et  al.,  2003],  These  ozone  mixing  ratios  are 
used  for  the  daytime  heating  and  cooling  rate  calculations. 

[17]  For  the  nighttime  cooling  rate  calculations,  these 
mixing  ratios  are  scaled  by  the  night-to-day  profile  in 
Figure  3b,  which  is  based  on  one-dimensional  photochem¬ 
ical  model  calculations  [Siskind  et  al.,  1995],  Use  of  a 
single  profile  reflects  the  fact  that  the  magnitude  of  the 
diurnal  ozone  variation  in  the  upper  stratosphere  and  lower 
mesosphere  does  not  show  much  seasonal  or  geographical 
variability  [e.g.,  Connor  et  al.,  1994;  Ricard  et  al.,  1996], 
This  correction  only  slightly  increases  our  nighttime  cooling 
rates,  since  the  cooling  rate  due  to  9.6  //m  ozone  emission  is 
<1  K  day~:  above  0.1  hPa  [Lopez-Puertas  and  Taylor, 
2001],  For  the  eclipse  simulations,  this  correction  was  not 
applied  within  umbral  shadow  regions. 

[is]  In  all  the  runs  reported  here,  C02  was  set  to  a  global 
constant  value  of  350  ppmv. 

3. 1.1.2.  Simplified  Model 

[19]  Since  CL1RAD  is  rigorously  validated  up  to  0.01  hPa 
only,  we  typically  run  NOGAPS-ALPHA  with  the  model 
top  at  0.005  hPa,  with  model  layers  at  0.01-0.005  hPa  a 
heavily  damped  sponge  region.  In  the  simulations  to  be 
described,  however,  we  extend  the  model  to  0.0005  hPa  in 
an  effort  to  study  possible  far-field  radiation  of  eclipse¬ 
generated  gravity  waves  to  higher  altitudes. 

[20]  The  dashed  curve  in  Figure  3c  shows  a  CLIRAD 
shortwave  heating  rate  calculated  off-line  at  a  specific 
location  in  the  eclipse  path  based  on  +8  hour  temperature 
and  constituent  profiles  from  a  NOGAPS-ALPHA  control 
run,  to  be  described  later.  At  altitudes  just  above  0. 1  hPa, 
the  reduction  with  height  of  the  ozone  heating  rates  abates, 
and  a  secondary  layer  of  radiative  heating  develops.  That 
heating  layer  extends  to  the  highest  model  altitudes. 

[21]  This  secondary  heating  rate  layer  is  problematic  for 
two  reasons.  First,  it  extends  to  altitudes  above  0.01  hPa 
(gray  curve  in  Figure  3  c)  where  CLIRAD  is  not  validated 
and  hence  the  heating  rates  are  unreliable.  Second,  our 
primary  focus  is  on  eclipse-induced  waves  generated  by 
reduced  shortwave  ozone  heating  in  the  stratosphere.  This 
secondary  high-altitude  layer  of  radiative  heating  can  pro¬ 
duce  additional  thermal  eclipse  forcing  that  can  potentially 
complicate  or  even  mask  the  wave  signals  radiating  upward 
from  the  reduced  ozone  heating  in  the  stratosphere. 

[22]  To  eliminate  this  secondary  heating  rate  layer  in  this 
study,  above  0.1  hPa  we  scaled  down  the  ozone  mixing 
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(a)  December  Ozone  Climatology  (ppmv)  (b)  Night-to-Day  Ratio  (c)  SW  Heating  Rate 


Figure  3.  (a)  Zonal-mean  ozone  mixing  ratios  (ppmv)  for  December  used  for  NOGAPS-ALPHA 
daytime  heating  and  cooling  rate  calculations,  (b)  Ozone  night-to-day  ratio,  used  to  scale  nighttime  ozone 
mixing  ratios  in  the  cooling  rate  calculations,  (c)  CL1RAD  radiative  heating  rates  at  indicated  location 
from  single  column  model  (dashed  curve).  Remaining  curves  show  changes  after  scaling  down  of  high- 
altitude  ozone  in  Figure  3a  using  equation  (1)  (dotted  curve)  and  then  further  removal  of  infrared  CO2 
heating  rates  (solid  curve).  Gray  line  shows  nominal  0.01  hPa  level,  below  which  CL1RAD  rates  are 
valid. 


ratios  in  Figure  3a  using  a  pressure-dependent  multiplicative 
coefficient 


/  logp  —  log  0.0 1  \2 

Viog0.001-iog0.01y  ’ 


(1) 


where  p  is  model  pressure  in  hPa,  and  e(p)  =  0  for  p  < 
0.001  hPa.  The  dotted  curve  shows  the  heating  rate  after 
imposing  this  ozone  reduction.  We  see  it  eliminates  nearly 
all  the  radiative  heating  at  p  <  0.01  hPa. 

[23]  From  0.1  to  0.01  hPa  a  residual  contribution  to  the 
secondary  heating  rate  layer  remains,  due  mostly  to  CCH 
absorption  in  the  near  infrared.  The  solid  curve  shows  the 
heating  rate  profile  after  CCH  heating  has  been  deactivated 
in  CLIRAD.  While  the  stratospheric  heating  is  reduced 
slightly,  the  secondary  layer  is  now  almost  completely 
eliminated. 

[24]  In  what  follows,  we  show  results  from  NOGAPS- 
ALPHA  runs  in  which  we  have  eliminated  CO2  heating  and 
reduced  the  high-altitude  O3  heating  in  this  way,  so  that  we 
can  focus  on  atmospheric  responses  to  reduced  stratospheric 
shortwave  heating  due  to  the  eclipse. 

3.1.2.  Hindcast  Runs 

[25]  Our  hindcast  experiments  use  a  T79L68  model 
formulation  extending  to  0.0005  hPa.  A  “cold  start”  ini¬ 
tialization  procedure  is  used  in  which  analyzed  winds, 
geopotential  heights  and  moisture  are  read  in  at  reference 
pressure  levels  on  a  1°  x  1°  grid  and  interpolated  to  the 
model’s  quadratic  Gaussian  grid  and  hybrid  a-p  levels. 


Initial  virtual  temperatures  are  computed  hydrostatically 
from  the  geopotentials.  The  model  was  then  forwarded  in 
time  without  assimilation  update  cycles  using  a  model  time 
step  of  300  s. 

[26]  From  1000  to  10  hPa,  these  initialization  fields  were 
specified  by  archived  NOGAPS  analysis  for  4  December 
2002  at  0000  UTC,  generated  operationally  at  the  Fleet 
Numerical  Meteorology  and  Oceanography  Center 
(FNMOC)  by  the  then-operational  Navy  multivariate  opti¬ 
mum  interpolation  (MVOI)  system  [ Barker ,  1992].  From  10 
to  0.4  hPa  we  use  FNMOC’s  “STRATOI”  analysis  of 
winds  and  geopotentials  [see  Goerrs  and  Phoebus,  1992, 
section  4].  At  altitudes  above  0.4  hPa,  where  there  are  no 
Navy  analysis  fields  for  this  date,  we  extrapolate  the  0.4  hPa 
STRATOI  fields  upward  by  progressively  relaxing  them 
with  increasing  altitude  to  zonal-mean  climatological  winds 
from  the  UARS  Reference  Atmosphere  Project  [Swinbank 
and  Ortland,  2003]  and  temperatures  from  the  1986 
COSPAR  International  Reference  Atmosphere  [. Fleming  et 
al.,  1990]:  for  algorithm  details,  see  Eckermann  et  ah 
[2004],  This  final  global  initial  state  within  NOGAPS- 
ALPHA  is  adjusted  for  hydrostatic  balance,  then  run 
through  a  nonlinear  normal  mode  initialization  procedure 
[Errico  et  ah,  1988]  to  improve  dynamical  balance  and 
suppress  spurious  gravity  wave  generation  due  to  unbal¬ 
anced  initial  conditions. 

[27]  At  the  surface,  ice  concentrations,  land/sea  surface 
temperatures  (SSTs)  and  snow  depths  are  initialized  using 
FNMOC  analysis.  Ice  and  SSTs  are  updated  from  archived 
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Figure  4.  (a)  Dispersion  curves  for  Lamb  waves  (thick  dot-dashed)  and  gravity  waves  at  their  turning 
points  ( m  =  0)  using  fully  compressible  (black  solid  curve)  and  pseudo-incompressible  (gray  solid  curve) 
dispersion  relations,  plotted  versus  horizontal  wavenumber  Ktot  and  intrinsic  frequency  Cj.  Horizontal 
dashed  lines  show  the  Cj  =  \f  \  and  Cj  =  N  limits,  (b)  Expansion  of  lower  left  comer  of  Figure  4a,  showing 
space-time  characteristic  of  eclipse  shadow  motion  for  V  =  700  m  s_1  (dotted  curve  with  diamonds)  and 
dispersion  curves  for  internal  gravity  waves  of  2n/\m\  =  100  km  (dot-dashed  curves  with  asterisks)  under 
both  fully  compressible  (black)  and  pseudo-incompressible  (gray)  dispersion  relations. 


analysis  every  24  hours:  grid  point  values  at  the  intervening 
times  are  linearly  interpolated. 

3.1.3.  Eclipse  Simulations 

[28]  The  hindcast  run  described  above  for  4  December 
2002  represents  our  control  simulation.  Our  eclipse  simula¬ 
tion  repeats  the  same  run,  keeping  everything  the  same  until 
045 1  UTC  when  the  eclipse  begins.  Then,  at  each  model  time 
step,  the  global  Euv{ A,  cp,  t)  field  at  the  current  universal  time 
t  of  the  model  simulation  is  read  in,  where  (A,  cp )  are 
longitude  and  latitude.  These  fields  scale  down  the  solar 
constant  used  in  the  radiative  heating  rate  calculations,  as 

S{\,<p,t)=S0[l-Euv(\,<P,t)].  (2) 

When  eclipsing  abates  at  1011  UTC,  the  simulation 
proceeds  as  before  with  Euv(\,  (p,  t)  =  0. 

[29]  By  simply  projecting  EUV(X,  cp,  t)  in  a  vertical 
column  through  the  atmosphere  via  equation  (2),  we  ignore 
vertical  variations  in  the  UV  obscuration  pattern.  These 
changes  with  height  will  be  largest  where  the  Sun  is  low 
in  the  sky  (cos  (j>  — >  0).  Since  Euv  values  in  Figure  1  are 
generally  small  near  the  terminator,  omission  of  these 
effects  here  should  not  introduce  significant  errors  given 
other  uncertainties  in  our  estimated  Euv{\,  (p,  t )  fields  (see 
Appendix  A). 

[30]  In  common  with  other  atmospheric  models, 
NOGAPS-ALPHA  usually  updates  its  radiative  heating 
and  cooling  rates  every  1-2  hours.  To  simulate  the  radiative 
response  of  the  atmosphere  to  rapidly  moving  eclipse 


shadows  via  equation  (2),  we  updated  the  radiative  heating 
and  cooling  rates  at  every  model  time  step,  both  in  the 
control  and  eclipse  simulations. 

[31]  In  both  runs,  instantaneous  global  model  fields  were 
saved  spectrally  at  one  hour  intervals. 

3.1.4.  Sensitivity  to  Anticipated  Dynamical  Responses 

[32]  Figure  4a  plots  horizontal  wavenumber  Ktot  versus 
intrinsic  frequency  Cj  within  the  nominal  internal  gravity 
wave  ranges  of  27r !Ktot  S;  10  km  and  f2  ^  Cj2  ^  N2,  where 
/  is  the  inertial  frequency  and  N  is  the  buoyancy  frequency. 
The  thick  solid  curve  shows  the  fully  compressible  non¬ 
hydrostatic  turning  point  curve,  where  vertical  wavenumber 
m  =  0  (Cj  =  Cjtp)  and  the  gravity  wave  reflects  vertically. 
Lamb  waves  occur  along  the  dot-dashed  curve,  which  also 
approximately  demarks  the  gravity  wave  turning  point  in 
the  hydrostatic  limit  [see  Marks  and  Eckermann,  1995, 
Figure  1],  These  two  curves  nearly  overlay  at  small  Ktot, 
then  begin  to  diverge  at  Ktot  S;  1-2  x  10-5  eye  m-1 
(wavelengths  ^50  km),  demarking  the  point  at  zonal  wave 
numbers  ^400  where  hydrostatic  dynamical  cores  such  as 
NOGAPS-ALPHA  start  exhibiting  errors  due  to  omitted 
nonhydrostatic  effects  [see,  e.g.,  Kasahara  and  Qian,  2000, 
Figures  2  and  4],  The  large-scale  eclipse  shadows  in  Figure  1 
are  expected  to  force  gravity  waves  with  horizontal  wave¬ 
lengths  of  thousands  of  kilometers  [Chimonas,  1970;  Fritts 
and  Luo,  1993],  placing  them  in  the  far  bottom  left  comer  of 
Figure  4a,  well  away  from  this  large  Ktot  region  where  the 
hydrostatic  approximation  fails.  This  justifies  our  use  of  a 
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Figure  5.  Fourier  solutions  of  the  gravity  wave  response  to  thermal  eclipse  forcing  (8)  based  on  the 
linearized  pseudo-incompressible  equations  of  Fritts  and  Luo  [1993].  First  column  plots  horizontal  cross 
sections  at  z  =  30  km  above  the  peak  in  eclipse  forcing  F  in  (a)  relative  potential  temperature,  (b)  zonal 
wind,  and  (c)  meridional  wind,  with  corresponding  vertical  cross  section  versus  x  along  the  eclipse 
shadow  axis  y  =  0  (gray  arrow)  for  (d)  relative  potential  temperature  and  (e)  zonal  wind,  (f)  Vertical  cross 
section  of  the  eclipse  heating  rate  perturbation  F,  normalized  by  its  peak  absolute  value  Q  =  —12  K 
day-1,  is  plotted.  Third  column  plots  horizontal  cross  sections  at  z  =  30  km  above  the  peak  F  of 
(g)  vertical  velocity  and  (h)  horizontal  divergence,  and  fourth  column  plots  vertical  cross  sections  along 
y  =  0  of  (i)  vertical  velocity  and  (j)  horizontal  divergence. 


T79  hydrostatic  global  spectral  model  to  simulate  eclipse- 
induced  wave  responses. 

[33]  Figure  4b  focuses  on  this  lower  left  portion  of 
Figure  4a.  Free-propagating  (internal)  gravity  waves  exist 
in  the  region  below  the  solid  black  turning  point  curve  (oj  < 
Cjtp )  and  above  the  dashed  curve  (<u  >  |/|).  The  dotted  line 
with  diamonds  shows  intrinsic  frequencies  for  a  wave 
disturbance  whose  horizontal  wavenumber  is  coaligned 
with  an  eclipse  shadow  velocity  V  =  700  m  s-1  and  whose 
horizontal  phase  speed  is  stationary  in  an  “eclipse  frame” 
moving  at  this  speed  V :  see  section  3.2  for  further  details. 
This  curve  also  represents  the  characteristic  space-time 
scales  of  the  thermal  stratospheric  forcing  due  to  the  eclipse. 
The  wave  field  it  generates  is  detennined  to  first  order  by 
how  directly  these  space-time  forcing  characteristics  project 
onto  various  wave  dispersion  curves  in  Figure  4b  [Salby 
and  Garcia,  1987].  We  see  that  this  forcing  curve  lies  above 
the  compressible  turning  point  and  Lamb  wave  curves 
(since  V  >  cs  ~  320  m  s-1,  where  cs  is  the  speed  of 
sound),  and  so  projects  most  directly  onto  a  vertically 
evanescent  (external)  gravity  wave  response.  The  wave 
dispersion  curves  closest  to  the  eclipse  forcing  curve  are 
the  (external)  Lamb  waves  and  high-frequency  (long  vertical 
wavelength)  internal  gravity  waves,  both  of  which  have  been 


predicted  by  theoretical  eclipse-forcing  models  [Chimonas, 
1970,  1973;  Fritts  and  Luo,  1993],  Vertically  deep  nonsta¬ 
tionary  gravity  wave  responses  are  also  predicted  by  simple 
thermal  forcing  models  where  space-time  scales  of  the 
source  and  gravity  waves  are  mismatched,  as  here  [e.g., 
Holton  et  al.,  2002]. 

[34]  A  global  hydrostatic  Eulerian  spectral  model  such  as 
NOGAPS-ALPHA  can  accurately  simulate  all  of  these 
anticipated  wave  responses  [see,  e.g.,  Davies  et  al.,  2003], 

3.2.  Fritts-Luo  Gravity  Wave  Model 

[35]  To  compare  our  NOGAPS-ALPHA  results  with 
existing  models  of  eclipse-induced  gravity  waves,  we  con¬ 
sider  the  equations  derived  and  solved  by  Fritts  and  Luo 
[1993].  They  begin  with  the  pseudo-incompressible  equa¬ 
tions  of  Durran  [1989],  which  filter  out  acoustic  and  Lamb 
waves,  then  simplify  further  using  the  hydrostatic  approx¬ 
imation,  since  nonhydrostatic  wave  responses  are  not  an¬ 
ticipated  (see  Figure  4).  Following  Chimonas  [1970],  Fritts 
and  Luo  [1993]  omit  rotation,  then  derive  linearized  per¬ 
turbation  forms  of  the  equations,  seeking  steady  (stationary) 
wave  solutions  in  a  coordinate  system  (x,  y,  z)  that  moves  in 
the  x-dircction  at  the  eclipse  shadow  speed  V,  which  is 
assumed  constant.  Background  winds  are  ignored,  although 
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(a)  Eclipse  Umbra:  0600  UTC  (b)  Eclipse  Umbra:  0700  UTC 


(c)  Eclipse  Umbra:  0800  UTC  (d)  Eclipse  Umbra:  0900  UTC 


cooling  rate  (K  day1)  cooling  rate  (K  day  ') 

Figure  6.  Gray  curves  show  difference  fields  in  total  radiative  heating  rate  between  eclipse  and  control 
simulations,  Q'(p),  at  (a)  0600  UTC,  (b)  0700  UTC,  (c)  0800  UTC  and  (d)  0900  UTC,  plotted  here  as 
(positive)  cooling  rates.  Solid  curves  show  corresponding  —QlcUPse(p)  profiles,  where  Q™upse  is  the 
shortwave  heating  rate  at  the  geographical  location  of  totality  in  the  absence  of  eclipse  obscuration 
effects.  These  latter  profiles  were  computed  off-line  with  a  single  column  model  using  profiles  of 
temperature,  water  vapor,  ozone,  albedo  and  ground  temperature  from  the  NOGAPS-ALPHA  control 
simulation  at  the  indicated  location  and  time.  Differences  lower  down  are  due  in  part  to  cloud  effects  not 
included  in  the  single  column  calculation. 


a  constant  flow  speed  could  be  added  to  V .  The  resulting 
equations  are: 


Vux  +  px/p  =  0, 

(3) 

Vvx  +Py/p  =  0, 

(4) 

o' 

II 

l<^> 

Qi 

to 

1 

(5) 

vy  +  wz-wg/c ;  =F/e, 

(6) 

V9X  +  6N2w/g  =  F, 

(7) 

where  u,  v,  w,  p  and  9  are,  respectively,  the  perturbations  of 
the  velocity  components,  pressure,  and  potential  tempera¬ 
ture.  The  mean  quantities  8  and  p  (atmospheric  density) 


are  functions  of  z  only.  The  speed  of  sound  is  cs,  the 
gravitational  acceleration  is  g,  and  the  mean  buoyancy 
frequency  is  N  =  (g9J9)112. 

[36]  For  the  eclipse-induced  diabatic  heating  rate  pertur¬ 
bation,  F,  Fritts  and  Luo  [1993]  chose  a  Gaussian  of  the 
form 


F  =  Qex  p 


(x2  +/)  z2  ' 

2al  2^1  ’ 


(8) 


with  oy=  1460  km,  cr,=  10  km,  and  \  Q\  =  12  K  day-1  (the  sign 
convention  for  Q  is  discussed  below).  Other  parameter  values 
in  the  Fritts-Luo  model  are:  N  =  0.02  s-1,  g  =  9.8  m  s"  2, 
T  =  2.7  x  10-2  km-1,  and  a  scale  height  of  7.8  km.  Note 
too  that  z  =  0  in  this  model  occurs  at  the  height  of  the  peak 
eclipse-induced  cooling  rate  \Q\,  not  at  the  surface. 

[37]  We  compute  numerical  solutions  in  the  t  — >  00  limit 
using  a  Fourier  method  outlined  in  Appendix  B.  The  first 
and  third  columns  of  plots  in  Figure  5  show  horizontal  cross 
sections,  at  z  =  30  km  above  the  peak  eclipse  forcing,  of 
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(a)  0.22  hPa 


(c)  8.67  hPa 


(b)  1.60  hPa 


(d)  46.84  hPa 


Figure  7.  Maps  of  radiative  cooling  rates  (K  day-1)  due  to  the  eclipse  at  0700  UTC  from  NOGAPS- 
ALPHA  runs,  plotted  at  (a)  0.22  hPa,  (b)  1.6  hPa,  (c)  8.7  hPa,  and  (d)  46.8  hPa.  Values  and  locations  of 
maxima  are  given  in  the  top  left  of  each  map. 


wave-induced  perturbations  of  relative  potential  temperature 
919,  zonal  wind  u,  meridional  wind  v,  vertical  velocity  w  and 
horizontal  divergence  D  =  ux  +  vy.  The  temperature  and 
velocity  plots  can  be  compared  to  Fritts  and  Luo  [1993, 
Figures  1  and  2],  which  span  the  much  smaller  horizontal 
subrange  of  —  1  x  104  km  <x  <  0  and  —0.6  x  104  km  <  y  < 
+0.6  x  104  km.  The  Fourier  solutions  within  this  ( x ,  y ) 
subrange  are  similar  in  form  and  magnitude  to  those  of 
Fritts  and  Luo  [1993],  apart  from  a  uniform  180°  phase 
difference.  This  sign  ambiguity  seems  to  originate  in  Q, 
which  is  set  to  12  K  day  1  by  Fritts  and  Luo  [1993], 
Flowever,  F  in  (6)  and  (7)  is  a  heating  rate  perturbation 
[Durran,  1989],  and  so  imposition  of  peak  eclipse-induced 


cooling  via  (8)  implies  Q  =  — 12  K  day-1,  a  sign  convention 
we  have  used  here. 

[38]  The  second  and  fourth  columns  of  plots  in  Figure  5 
show  corresponding  vertical  cross  sections  of  the  wave 
fields  along  the  y  =  0  eclipse  shadow  axis.  The  exception 
is  Figure  5f,  which  profdes  the  normalized  eclipse-induced 
cooling  rate  function  FIQ,  since  the  v  field  beside  it  is 
antisymmetric  about  y  =  0  and  thus  its  vertical  cross  section 
is  featureless.  These  cross  sections  all  show  a  vertically 
deep  wave  response  (vertical  wavelengths  ~  100  -200  km), 
which  grows  in  amplitude  with  increasing  altitude. 

[39]  Our  Fourier  solutions  showed  considerable  sensitiv¬ 
ity  to  our  choice  for  T.  The  solutions  in  Figure  5  were 
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Figure  8.  Eclipse-induced  radiative  cooling  rates  from 
NOGAPS-ALPHA  runs  at  0700  UTC,  averaged  over  the 
indicated  40°  longitude  and  20°  latitude  band  centered 
about  the  eclipse  umbra  location  at  this  time. 

derived  using  F  =  2.45  x  10  2  km-1,  a  value  about  10% 
smaller  than  that  used  by  Fritts  and  Luo  [1993],  which  gave 
better  agreement  with  their  results.  The  reason  for  this 
sensitivity  can  be  gleaned  from  Figure  4b,  whose  gray  solid 
curve  shows  the  gravity  wave  turning  point  curve  under  the 
pseudo-incompressible  dispersion  relation  (B6)  with  rota¬ 
tion  retained  ( m  =  0,  u>  —  CjPiTP).  We  see  that  the  pseudo- 
incompressible  turning  point  frequencies  Cjp1tp  exceed  the 
fully  compressible  values  uTP  (black  solid  curve),  thus 
permitting  a  greater  range  of  harmonic  pairs  ( Ktot ,  CS)  to 
be  freely  propagating  gravity  waves.  Furthermore,  the 
space-time  characteristic  curve  of  eclipse  motion  (dotted 
line  with  diamonds)  essentially  overlays  this  pseudo-incom¬ 
pressible  turning  point  curve.  Thus  slight  changes  in  V  or  F 
which  change  the  relative  positions  of  these  two  curves 
allow  for  significant  increases  or  decreases  in  forcing  of 
freely  propagating  gravity  wave  harmonics  in  this  model. 
More  importantly,  Figure  4b  shows  that  the  pseudo-incom¬ 
pressible  solutions  will  overestimate  internal  gravity  wave 
generation  by  the  moving  eclipse  shadow,  by  allowing 
supersonic  wave  phase  speeds  to  propagate  vertically  as 
internal  gravity  waves.  This  reflects  a  breakdown  in  the 
pseudo-incompressible  approximation,  which  requires  the 
Lagrangian  timescales  of  these  forced  wave  disturbances  to 
be  much  longer  than  timescales  for  sound  wave  propagation 
[Durran,  1989],  or  equivalently,  horizontal  phase  speeds  to 
be  much  slower  than  cs.  As  Figure  4b  attests,  this  criterion  is 
not  satisfied  generally  for  thermal  forcing  by  supersonically 
moving  eclipse  shadows. 

[40]  Despite  these  caveats,  these  pseudo-incompressible 
Fritts  and  Luo  [1993]  solutions  are  valuable  in  providing 
guidance  on  the  three-dimensional  shapes  and  amplitudes  of 
stationary  gravity  wave  responses  anticipated  from  broadly 


realistic  eclipse-induced  thennal  cooling  of  the  stratosphere, 
for  cross  comparison  with  the  NOGAPS-ALPHA  results. 

4.  NOGAPS-ALPHA  Results 

4.1.  Eclipse  Radiative  Cooling  Rate  Footprint 

[41]  In  any  given  NOGAPS-ALPHA  grid  box,  the  net 
radiative  heating  rate  profile  Q(p)  =  Qsw{p )  +  Qlw(p),  where 
Qs"  is  the  shortwave  solar  heating,  Q"  is  the  longwave 
cooling  contribution  and  p  is  model  pressure.  Gray  curves  in 
Figure  6  plot  vertical  profiles  of  Q'  =  QecUpse  -  QC0nt,0i,  the 
difference  in  these  net  heating  rate  profiles  between  eclipse 
and  control  simulations,  plotted  at  umbral  shadow  locations 
from  0600  UTC  to  0900  UTC  (see  Figure  la).  These 
profiles  all  reveal  uniformly  negative  heating  rate  perturba¬ 
tions  Q'  due  to  the  eclipse  UV  obscuration  E(m  and  so  are 
plotted  in  Figure  6  as  (positive)  cooling  rates. 

[42]  Within  these  regions  of  totality,  Ejjv  =  1  and  thus 
from  (2)  shortwave  heating  is  totally  shut  off  (Q?<£upse  =  0). 
Thus,  if  the  longwave  cooling  does  not  change  significantly 
between  eclipse  and  noneclipse  conditions,  then  Q '  k, 
— (/control-  To  test  this,  the  black  curves  in  Figure  6  plot 
off-line  estimates  of  —  Control  using  saved  output  profiles 
from  the  control  simulation  at  these  times  and  locations. 
These  —  Q^mroi  profiles  essentially  overlay  the  stratospheric 
Qr  profiles. 

[43]  The  gray  £/  profiles  in  Figure  6  occur  within  totally 
eclipsed  atmospheric  regions.  In  penumbral  regions  (0  < 
Euv<  1),  reduced  (rather  than  zero)  solar  shortwave  heating 
occurs  during  eclipse  passages,  and  so  the  radiative  cooling 
rate  signature  of  the  eclipse  is  reduced  here.  To  define  the 
overall  horizontal  structure  of  these  eclipse-induced  cooling 
rate  footprints,  Figure  7  plots  maps  of  Q'  from  the 
NOGAPS-ALPHA  runs  at  0700  UTC  at  4  different  pressure 
levels,  ranging  from  the  lower  mesosphere  (0.22  hPa)  to  the 
lower  stratosphere  (47  hPa).  The  horizontal  structures 
generally  resemble  that  of  the  corresponding  0700  UTC 
Ei jy  footprint  in  Figure  If. 

[44]  Maps  of  tropospheric  radiative  heating  rate  perturba¬ 
tions  like  those  in  Figure  7  (not  shown)  reveal  much  greater 
spatial  inhomogeneity,  apparently  related  to  cloud  fields 
simulated  in  the  model.  To  assess  the  potential  for  a  net 
large-scale  forcing  effect  [e.g.,  Chimonas,  1973],  Figure  8  plots 
a  mean  0700  UTC  Q'(p)  profile  averaged  over  a  broad 
geographical  region  centered  about  the  umbra.  It  reveals 
smaller  but  significant  eclipse-induced  radiative  cooling 
rates  throughout  the  troposphere,  with  forcing  peaking  in 
this  case  near  850  hPa. 

4.2.  Middle  Atmospheric  Responses 

[45]  We  focus  initially  on  the  atmospheric  responses  at 
0800  UTC,  since  they  are  morphologically  similar  (but  of 
larger  amplitude  and  size)  to  those  simulated  at  earlier 
eclipse  times.  By  0800  UTC  the  eclipse  shadow  is  moving 
near  its  limiting  speed  of  ^-700  m  s  1  over  oceanic  regions 
where  the  Sun  is  high  in  the  sky  (Figures  lh  and  2).  Figure  9 
plots  maps  of  difference  fields  between  eclipse  and  control 
runs  at  1.6  hPa,  a  level  near  the  peak  in  the  eclipse-induced 
radiative  cooling  rate,  and  at  0.01  hPa,  a  region  ~30— 40  km 
above  the  peak  eclipse  cooling  in  a  region  where  no  local 
thermal  forcing  due  to  the  eclipse  is  anticipated  in  these  runs 
(see  Figure  6). 
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Figure  9.  Polar  orthographic  maps  of  modeled  atmospheric  responses  to  eclipse  passage  at  0800  UTC 
at  0.01  hPa  (columns  1  and  3)  and  1.6  hPa  (columns  2  and  4).  Parameters  plotted  are  (a-d)  temperature 
T'  and  vertical  velocity  W',  (e-h)  zonal  wind  If  and  horizontal  divergence  O',  and  (i-1)  meridional  wind 
V'and  geopotential  height  Z'.  The  minimum  and  maximum  values,  contour  interval  and  units  are  given  at 
the  bottom  of  each  map.  Positive  values  have  black  contours,  and  negative  values  have  thicker  gray 
contours.  Dashed  line  shows  passage  of  umbral  shadow,  extrapolated  to  longitudes  west  of  Africa.  Solid 
black  circle  shows  0800  UTC  location  of  the  umbra.  15x7  point  longitude-latitude  smoothing  was 
applied  to  these  fields. 


[46]  The  horizontal  structures  of  the  temperature  responses, 
T’  =  Teclipse  —  Tcontroi,  differ  notably  at  each  altitude.  At 
1.6  hPa  the  response  is  mostly  a  uniform  decrease  that  peaks 
at  T'  =  —1.0  K  slightly  behind  the  umbral  shadow.  This 
reflects  the  direct  atmospheric  cooling  effect  of  reduced 
shortwave  ozone  heating  of  this  model  layer  due  to  obscu¬ 
ration  of  the  solar  disc.  The  integrated  effect  of  this 
atmospheric  temperature  decrease  is  also  seen  in  the  geo¬ 
potential  heights,  7f  which  show  a  similar  structure  to  the 
1.6  hPa  T'  field. 

[47]  The  “far  field”  temperature  response  at  0.01  hPa  is 
completely  different.  It  is  characterized  by  a  warm  central 
core,  peaking  at  0.22  K  over  Africa,  surrounded  by  a  belt  of 
weaker  negative  T'  values.  Similar  annular  responses  are 
simulated  at  both  altitudes  in  vertical  velocities,  W',  and 
horizontal  divergences,  D' . 


[48]  The  core  of  negative  D'  and  W'  values  surrounded  by 
a  belt  of  weaker  positive  values  implies  a  convergent 
circulation  in  which  air  moves  inward  toward  a  central 
point  located  near  southern  Africa,  then  descends,  with 
weaker  divergence  and  ascent  at  the  outer  edges.  This  is 
reflected  in  a  zonal  velocity  response  U'  that  is  roughly 
zonally  antisymmetric  about  the  20°E  meridian  near  Africa: 
eastward  of  this  meridian,  a  lobe  of  negative  U'  values 
occurs,  peaking  near  the  umbral  shadow,  while  to  the  west 
of  this  meridian  a  lobe  of  positive  U'  values  occurs  well 
away  from  the  moving  eclipse  shadow.  Similarly,  the 
meridional  velocity  response  V'  is  approximately  antisym¬ 
metric  about  a  latitude  circle:  poleward,  a  lobe  of  positive  V' 
values  occurs,  and  equatorward,  a  lobe  of  negative  V  values 
occurs.  The  peak  magnitudes  of  all  four  horizontal  velocity 
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(c)  Vertical  Velocity  0800  UTC  (d)  Geopotential  Height  0800  UTC 


Figure  10.  Vertical  cross  sections  of  0800  UTC  modeled  atmospheric  responses  to  eclipse  passage 
along  umbral  trajectory  (dashed  curve  in  1.6  hPa  plots  of  Figure  9):  (a)  temperatures,  contour  interval 
0.2  K,  with  ±0.02  K,  ±0.05  K  and  ±0. 1  K  contours  also  shown;  (b)  zonal  wind,  contour  interval  0.2  ms-1; 
(c)  vertical  velocity,  contour  interval  0.2  cm  s-1;  and  (d)  geopotential  height,  contour  interval  10  m.  In  all 
plots,  positive  values  have  black  contours,  and  negative  values  have  thicker  gray  contours.  Horizontal 
dotted  lines  mark  0.01  hPa  and  1.6  hPa  levels  mapped  in  Figure  9.  Vertical  dashed  line  shows  0800  UTC 
umbral  shadow  location. 


lobes  are  about  the  same,  roughly  0.5  ms-1  at  1.6  hPa  and 
1  m  s_1  at  0.01  hPa. 

[49]  To  study  the  vertical  structure  of  these  responses, 
Figure  10  plots  altitude  cross  sections  along  the  extrapolated 
umbral  path  depicted  by  the  dashed  curves  in  the  1.6  hPa 
plots  of  Figure  9.  Temperature  cross  sections  show  a  broad 
region  of  cold  upper  stratospheric  anomalies  due  to  eclipse- 
induced  radiative  cooling  of  this  region  (Figure  6).  At 


altitudes  above  0.1  hPa,  the  response  transitions  to  a  weak 
warming  in  the  wake  of  the  eclipse  shadow,  giving  way  to  a 
weaker  cooling  response  near  the  Greenwich  meridian  that 
tilts  westward  with  height.  On  descending  through  the 
stratosphere,  the  mean  cooling  response  progressively 
weakens  and  disappears,  giving  way  to  a  weak  warming 
peak  near  70  hPa  before  reappearance  of  temperature 
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Temperature  0900  UTC 


0900  UTC 


Vertical  Velocity  0900  UTC 


Figure  11.  Same  presentation  as  Figure  9  but  plotting  the  modeled  atmospheric  response  to  the  eclipse 
at  0900  UTC. 


decreases  in  the  troposphere  due  to  direct  radiative  cooling 
of  these  layers  (Figure  8). 

[so]  The  two-lobed  zonal  wind  response  in  Figure  10  is 
vertically  very  deep.  The  positive  {If  >  0)  lobe  tilts 
noticeably  to  the  west.  The  negative  lobe  tilts  slightly 
eastward  with  height,  though  the  if  responses  at  and  ahead 
of  the  umbral  shadow  (vertical  dashed  curve)  are  almost 
columnar  ( U l  =  0),  as  are  the  structures  in  vertical  velocities 
and  geopotential  heights  at  these  locations.  Apart  from  the 
increase  in  overall  size  and  amplitude,  the  most  notable 
change  in  0800  UTC  fields  compared  to  earlier  times  is  greater 
eastward  tilting  with  height  of  the  middle-atmospheric 
responses.  For  example,  the  vertical  cross  sections  of  7f  in 
Figure  lOd  reveal  greater  eastward  tilting  with  height  of  the 
response  behind  the  umbral  shadow,  associated  with  forma¬ 
tion  of  an  isolated  7'  minimum  at  high  altitudes  just  east 
of  the  Greenwich  meridian.  The  horizontal  0.01  hPa  Z' 
map  in  Figure  9k  shows  this  is  due  to  formation  of  a  weak 
V-shaped  response. 

[51]  By  0900  UTC,  bow-wave-like  responses  become 
much  stronger.  For  example,  Figure  11  reveals  7'  responses 


at  both  altitudes  with  much  clearer  V-shaped  structure, 
and  formation  of  weak  positive  anomalies  over  Africa  at 
0.01  hPa.  At  that  same  location  and  altitude,  positive  W' 
and  D'  anomalies  also  arise.  These  bow-wave  responses 
distort  the  previously  symmetric  zonal  winds  responses  if, 
particularly  at  higher  altitudes  (Figure  Ilf).  Associated  with 
this  is  greater  eastward  tilting  of  the  vertical  cross  sections 
of  U',  W'  and  7'  in  Figure  12.  The  direct  thermal  cooling 
signatures  of  the  eclipse  in  the  stratosphere  and  troposphere 
in  Figure  12a  are  now  separated  by  a  distinct  region  of  weak 
positive  temperature  anomalies  in  the  lower  stratosphere 
and  mesosphere  behind  the  eclipse  shadow,  indicating  both 
upward  and  downward  propagation  of  bow-wave-like 
responses  from  upper  stratospheric  thermal  eclipse  forcing. 

[52]  Small-scale  horizontal  structure  is  seen  in  W'  in 
Figure  12c  over  the  African  longitudes.  This  seems  to 
originate  from  changes  in  the  tropospheric  shortwave  cloud 
forcing  over  the  African  land  mass,  which  modifies  the 
convective  clouds  and  the  properties  of  explicitly  resolved 
convectively  generated  short-wavelength  gravity  waves, 
which  propagate  to  high  altitudes.  Their  effects  are  evidenced 


13  of  22 


D14105 


ECKERMANN  ET  AL.:  MODELING  ATMOSPHERIC  EFFECTS  OF  A  SOLAR  ECLIPSE 


D14105 


(a)  Temperature  0900  UTC 


(b)  Zonal  Wind  0900  UTC 


0  50  100 

(c)  Vertical  Velocity  0900  UTC 


0  50  100 


(d)  Geopotential  Height  0900  UTC 


Figure  12.  Same  presentation  as  Figure  10  for  the  0900  UTC  modeled  atmospheric  response  to  eclipse 
passage:  (a)  temperatures,  contour  interval  0.2  K,  with  ±0.02  K,  ±0.05  K  and  ±0.1  K  contours  also 
shown;  (b)  zonal  wind,  contour  interval  0.2  ms-1;  (c)  vertical  velocity,  contour  interval  0.2  cm  s-1;  and 
(d)  geopotential  height,  contour  interval  10  m.  Vertical  dashed  line  shows  0900  UTC  umbral  shadow 
location. 


by  the  noisier  nature  of  the  W'  and  D'  fields  in  Figure  11, 
even  after  significant  horizontal  smoothing. 

[53]  As  the  eclipse  nears  its  end  by  1000  UTC  (Figure  11), 
Figure  13  reveals  bow  wave-like  responses  dominating  at 
high  altitudes.  For  example,  the  symmetric  two-lobe  zonal 
wind  response  at  earlier  times  has  been  largely  replaced  in 
Figure  13b  by  a  well-defined  V-shaped  pattern  that  is 
approximately  symmetric  about  the  umbral  shadow  path 
(dotted  green  line).  Similar  V-shaped  responses  are  also 
seen  in  temperature,  vertical  velocity  and  geopotential 
height  in  Figures  13a,  13c,  and  13d. 


4.3.  Tropospheric  Responses 

[54]  Figure  12a  reveals  midtropospheric  cooling  of 
~0.1  K  due  to  local  radiative  cooling  within  the  eclipse 
shadow  (Figure  8).  Near  the  surface,  however,  much  larger 
atmospheric  temperature  drops  have  been  measured  over 
land  during  eclipses  due  to  radiative  surface  cooling  and 
accompanying  changes  in  turbulent  fluxes  of  latent  and 
sensible  heat  [e.g.,  Segal  et  al.,  1996].  To  investigate  this  in 
NOGAPS-ALPFIA,  Figure  14  plots  the  air  temperature 
response  T's  at  0700  UTC  in  the  lowest  full  model  layer. 
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(a)  Temperature,  0.01018  hPa,  1000  UTC  (b)  Zonal  Wind,  0.01018  hPa,  1000  UTC 


Figure  13.  Maps  of  modeled  atmospheric  responses  to  eclipse  at  1000  UTC  at  0.01  hPa  in  (a)  temperature, 
(b)  zonal  wind,  (c)  vertical  velocity,  and  (d)  geopotential  height.  Scales  and  units  are  given  in  the  color 
bar,  and  minimum  and  maximum  values  at  the  bottom  of  each  plot.  Positive  (negative)  values  have  white 
solid  (dotted)  contours.  These  contours  are  every  0.1  K  with  ±0.05  K  contours  also  shown  (Figure  13a), 
±0.2ms"1,±0.4ms“1,±l  ms~'  and±1.5  ms-1  (Figure  13b), ±0.25  cms-1  and ±0.5  cms-1  (Figure  13c), 
and  every  20  m  with  ±10  m  contours  also  shown  (Figure  13d).  Thick  dashed  line  shows  passage  of 
umbral  shadow,  and  solid  circle  shows  its  final  location  just  after  0900  UTC  (see  Figure  1).  15  x  7  point 
longitude-latitude  smoothing  was  applied  to  these  fields. 


Surface  Pressure  1000  UTC 


Surface  Temperature  0700  UTC 


Figure  14.  Modeled  surface  air  temperature  response  T's  at 
0700  UTC.  Scales  and  units  are  given  in  the  color  bar,  and 
minimum  and  maximum  values  are  given  at  the  bottom  of 
the  map.  Thick  dashed  line  shows  passage  of  umbral 
shadow,  and  solid  circle  shows  its  current  location. 


Figure  15.  Modeled  surface  pressure  response  p's  at 
1000  UTC.  Scales  and  units  are  given  in  the  color  bar,  and 
minimum  and  maximum  values  are  given  at  the  bottom  of 
the  map.  Positive  (negative)  values  have  white  solid  (dotted) 
contours:  values  are  ±0.05  hPa,  ±0.1  hPa,  0.3  hPa  and 
0.5  hPa.  Thick  dashed  line  shows  passage  of  umbral 
shadow,  and  solid  circle  shows  its  final  location  just  after 
0900  UTC  (see  Figure  1).  15  x  7  point  longitude-latitude 
smoothing  was  applied  to  these  fields. 
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Table  1.  Differences  Between  Fritts  and  Luo  [1993]  and 
NOGAPS-ALPHA  Wave  Solutions  in  Response  to  Thermal 
Stratospheric  Eclipse  Forcing 


Parameter 

Fritts  and  Luo  [1993] 

NOGAPS-ALPHA 

Dispersion  relation 

Major  Differences 
pseudo-incompressible 

compressible 

Shadow  speed  V 

constant 

time-varying3 

Eclipse  cooling 

constant*5 

variable0 

Peak  cooling 

12  K  day-' 

~27  K  day-1 

Gravity  wave  solutions 

stationary  in  V  frame 

no  restriction 

Lamb  waves 

no 

yes 

Times 

t  — >  oo  limit 

hourly 

Wind/shear 

Minor  Differences 
no 

yes 

Geometry 

Cartesian 

spherical 

Rotation 

no 

yes 

aSee  Figure  2. 
bSee  equation  (8). 
°See  Figures  1  and  6. 


located  just  13  m  above  the  ground  for  surface  pressures  of 
~  1 000  hPa.  We  see  atmospheric  cooling  of  up  to  ~4  K  that 
is  confined  over  land  (southern  Africa  and  Madagascar). 
These  values  are  within  the  2-10  K  range  reported  in  both 
near-surface  measurements  [e.g.,  Anderson  et  al.,  1972; 
Segal  et  al .,  1996;  Eaton  et  al.,  1997]  and  mesoscale  model 
simulations  [ Gross  and  Hense,  1999;  Vogel  et  al.,  2001] 
during  previous  eclipse  passages. 

[55]  Figure  15  plots  the  1000  UTC  surface  pressure 
response  p's.  Unlike  the  corresponding  T's  map  (not  shown), 
p's  exhibits  a  coherent  large-scale  response  over  both  land 
and  ocean  that  is  very  similar  in  overall  form  to  the  bow- 
wave  structures  encountered  at  high  altitudes  in  Figure  13  at 
this  time.  Amplitudes  here  are  ^0.1 -0.5  hPa,  having 
increased  with  time  during  the  eclipse  passage,  with  the 
peak  values  clustered  over  the  southern  African  land  mass. 

5.  Discussion 

[56]  We  now  compare  our  simulated  atmospheric  responses 
to  this  solar  eclipse  in  NOGAPS-ALPHA  with  some  rele¬ 
vant  observations  and  theories. 

5.1.  Stratospheric  Bow  Wave  Theory 

[57]  The  primary  goal  of  the  NOGAPS-ALPHA  runs  was 
to  see  whether  reduced  stratospheric  ozone  heating  during 
the  eclipse  generated  a  large-scale  internal  gravity  wave  at 
higher  altitudes  with  the  bow-wave  form  and  amplitudes 
predicted  by  the  linear  models  of  Chimonas  [1970]  and 
Fritts  and  Luo  [1993],  To  facilitate  this  comparison,  we 
regenerated  the  Fritts  and  Luo  [1993]  gravity  wave  solu¬ 
tions  in  Figure  5.  Exact  correspondences  between  these 
solutions  and  our  NOGAPS-ALPHA  results  are  not 
expected  because  of  a  number  of  differences  between  the 
two  formulations,  the  more  important  of  which  are  summa¬ 
rized  in  Table  1. 

[58]  For  this  particular  eclipse,  for  example,  the  NOGAPS- 
ALPFIA  runs  gave  peak  eclipse-induced  radiative  cooling 
rates  of  ~27  K  day-1 :  see,  e.g.,  Figures  6c  and  7b.  The  bow- 
wave  model  of  Fritts  and  Luo  [1993]  adopted  a  constant  peak 
cooling  rate  \Q\  of  12  K  day’  1  (see  section  3.2),  while 
Chimonas  [1970]  used  ~9  K  day-1  [see  Murgatroyd  and 


Goody,  1958;  Davies,  1982],  Thus  our  NOGAPS-ALPHA 
simulations  for  this  one  specific  eclipse  indicate  that  earlier 
model  studies  have  underestimated  the  strength  of  this 
peak  eclipse-induced  stratospheric  thermal  forcing  by  a 
factor  of  2-3. 

[59]  Atmospheric  responses  in  NOGAPS-ALPHA  appeared 
at  high  altitudes  quickly  and  rapidly  expanded  both  merid¬ 
ionally  and  zonally.  This  indicates  additional  forcing  of 
gravity  waves  that  are  nonstationary  with  respect  to  eclipse 
motion  and  thus  not  captured  by  the  Fritts  and  Luo  [1993] 
model.  These  waves  arise  because  of  more  realistic  space- 
time  variations  in  thermal  eclipse  forcing  (e.g.,  variable  V 
and  Q)  in  the  NOGAPS-ALPHA  runs  [e.g.,  Holton  et  al., 
2002],  As  a  result,  the  simulated  NOGAPS-ALPHA 
responses  at  early  times  showed  near-columnar  annular  or 
symmetric  responses  (Figures  9  and  10)  that  differed  from 
bow-wave  predictions.  However,  by  1000  UTC  these  early  - 
stage  responses  transitioned  to  new  forms  more  reminiscent 
of  bow-wave  model  predictions.  This  is  consistent  with  the 
long  time  limit  ( t  — >  00)  of  the  Fritts  and  Luo  [1993] 
solutions:  fully  compressible  vertical  group  velocity  calcu¬ 
lations  based  on  wavelength-frequency  pairs  in  Figure  4b 
indicate  that  even  the  very  long  vertical  wavelength  internal 
gravity  wave  modes  can  take  several  hours  to  propagate  to 
altitudes  30-40  km  above  the  stratospheric  eclipse  forcing. 

[60]  At  these  later  times  there  are  clear  similarities  between 
the  0.01  hPa  NOGAPS-ALPHA  responses  and  the  Fritts  and 
Luo  [1993]  solutions  in  Figure  5.  For  example,  Figure  5c 
predicts  an  antisymmetric  meridional  wind  response  V' 
about  the  y  =  0  axis  of  the  eclipse  shadow  motion,  distorted 
downstream  by  bow-wave  effects,  with  negative  (positive) 
perturbations  to  the  north  (south)  and  peak  magnitudes  of 
~0.25  ms-1  at  30  km  above  the  eclipse  forcing.  The 
NOGAPS-ALPHA  V'  response  at  0900  UTC  (Figure  lli) 
has  a  similarly  phased  antisymmetric  form  about  the  umbral 
shadow  trajectory,  but  has  larger  peak  magnitudes  of 
~1.7  m  s  1 .  The  Fritts  and  Luo  [1993]  zonal  velocity 
response  if  in  Figure  5b  is  symmetric  about  y  =  0,  with  a 
strong  negative  lobe  ahead  of  the  eclipse,  and  a  weaker 
V-shaped  positive  lobe  behind  it.  NOGAPS-ALPHA 
1000  UTC  fields  in  Figure  13b  show  a  similar  structure, 
though  with  less  difference  between  the  leading  and  lagging 
lobe  amplitudes.  The  Fritts  and  Luo  [1993]  temperature 
response  in  Figure  5a  is  symmetric  about  y  =  0,  dominated 
by  a  warm  anomaly  in  the  wake  of  the  eclipse,  peaking  at  a 
relative  amplitude  of  about  10-3:  for  a  mean  0.01  hPa 
temperature  of  ~200  K,  this  yields  T'  ~  0.2  K.  Warm 
symmetric  V-shaped  anomalies  in  the  wake  of  the  eclipse 
are  also  seen  in  NOGAPS-ALPHA  fields,  with  a  peak 
magnitude  of  ~0. 3-0.4  K  (Figures  11a  and  13a). 

[61]  The  vertical  velocity  and  divergence  responses  in 
Figures  5g  and  5h,  respectively,  are  predicted  to  have  a 
symmetric  three-lobe  structure,  with  weak  positive  anoma¬ 
lies  ahead  of  the  eclipse,  a  strong  V-shaped  negative 
anomaly  immediately  behind  the  umbra  and  a  third  weaker 
V-shaped  positive  anomaly  forming  further  behind  the 
umbra  (see  also  Figure  5j).  These  two  symmetric  lagging 
lobes  are  both  reproduced  in  the  NOGAPS-ALPHA  W'  and 
D '  fields.  For  example.  Figures  11c  and  1 1  g  show  this 
second  positive  lagging  lobe  having  just  formed  over  Africa 
within  the  main  negative  lagging  lobe.  Both  positive  and 
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negative  anomalies  become  larger  and  more  V-shaped  by 
1000  UTC  (Figure  13c). 

[62]  The  vertical  cross  sections  of  the  Fritts  and  Luo 
[1993]  solutions  in  Figure  5  all  show  vertically  deep 
responses  that  tilt  eastward  with  increasing  height,  consis¬ 
tent  with  long  vertical  wavelength  gravity  waves  that  are 
stationary  in  a  frame  moving  eastward  at  speed  V.  Deep 
eastward  tilted  responses  are  also  seen  in  the  0900  UTC 
NOGAPS-ALPF1A  velocity  and  geopotential  height  cross 
sections  in  Figure  12.  Similar  structure  is  masked  in  the 
temperature  cross  sections  in  Figure  12a  because  of  a  larger 
mean  temperature  decrease  in  the  upper  stratosphere  due  to 
local  radiative  cooling. 

[63]  In  summary,  the  mesospheric  NOGAPS-ALPFIA 
fields  at  0900-1000  UTC  show  clear  evidence  of  a  large- 
scale  V-shaped  internal  gravity  wave  forced  by  radiative 
cooling  of  the  stratosphere  within  the  eclipse  shadow.  The 
overall  horizontal,  vertical  and  amplitude  structures  agree 
broadly  with  the  predictions  of  Fritts  and  Luo  [1993].  In 
particular,  our  results  confirm  their  findings  that  mesospher¬ 
ic  wave  amplitudes  are  weak,  which  may  explain  the 
difficulty  in  observing  this  wave  response  during  eclipse 
passages. 

5.2.  Middle  Atmospheric  Observations 

[64]  Some  rocket  soundings  of  the  middle  atmosphere 
have  reported  temperature  decreases  in  the  5-12  K  range  at 
50-60  km  altitude  during  eclipse  passages  [ Ballard  et  al., 
1969;  Quiroz  and  Henry,  1973;  Randhawa,  1974;  Schmidlin 
and  Olsen,  1984].  Quiroz  and  Henry  [1973]  and  Schmidlin 
and  Olsen  [1984]  also  reported  substantial  increases  in 
meridional  wind  speeds,  peaking  at  20-40  m  s-1  at 
~60  km,  which  Quiroz  and  Henry  [1973]  interpreted  as  a 
balanced  circulation  response  to  eclipse-induced  changes  in 
the  lateral  temperature  gradients.  Our  NOGAPS-ALPHA 
simulations  were  unable  to  reproduce  these  observations. 
We  found  no  evidence  of  wind  or  temperature  responses 
anywhere  near  this  large  in  our  model  fields  at  50-60  km: 
temperature  decreases  here  were  ^1  K,  and  horizontal  wind 
changes  were  a  few  meters  per  second  at  most  (see,  e.g., 
Figure  12). 

[65]  Several  points  are  worth  noting.  First,  rocketsondes 
measure  instantaneous  profiles  with  small-scale  gravity  wave 
perturbations  superimposed,  and  the  typical  r.m.s.  ampli¬ 
tudes  of  these  oscillations  at  50-60  km  are  ~  1 0  m  s-1  for 
horizontal  winds  and  ^2-5  K  for  temperature  [Eckermann 
et  al.,  1995].  Given  wave  periods  as  short  as  5-10  min, 
serendipitous  phasing  among  gravity  waves  in  sequential 
rocket  profiles  during  eclipse  passages  could  yield  wind 
and  temperature  changes  of  the  order  of  those  reported,  but 
which  would  have  no  physical  connection  to  the  eclipse. 
Since  this  phasing  argument  should  just  as  easily  produce 
warming,  it  cannot  explain  why  all  the  studies  have 
observed  strong  cooling,  although  there  are  other  obser¬ 
vations  that  have  reported  no  measurable  changes  in 
middle  atmospheric  winds  or  temperatures  during  eclipses 
[Randhawa,  1973;  Ball  et  al.,  1980]. 

[66]  More  detailed  observations  would  clearly  help.  For 
example,  ground-based  temperature  and  wind  (Doppler) 
lidars  with  day-night  capabilities  [e.g.,  Alpers  et  al.,  2004] 
could  acquire  data  with  the  necessary  time-height  resolution 


to  separate  mean  and  gravity  wave  components  at  these 
altitudes  during  a  solar  eclipse. 

5.3.  Bow  Wave  Response  in  Surface  Pressure 

[67]  The  NOGAPS-ALPHA  surface  pressure  response  p's 
at  1000  UTC  in  Figure  15  exhibits  a  large-scale  bow-wave- 
like  structure  very  similar  to  that  seen  in  the  high-altitude 
responses  in  Figure  13.  Stratospheric  eclipse  cooling  should 
radiate  gravity  waves  to  both  higher  and  lower  altitudes 
[Chimonas,  1970],  For  the  March  1970  eclipse,  Chimonas 
and  Hines  [1970]  predicted  such  downward  radiation  of 
gravity  waves  should  yield  surface  pressure  oscillations 
p's  ~  ±0.01  hPa.  These  values  are  an  order  of  magnitude  or 
more  smaller  than  our  simulated  amplitudes  in  Figure  15, 
which  are  in  the  0.1 -0.5  hPa  range.  However,  these 
amplitudes  are  within  the  ±0.1-1  hPa  range  observed  in 
microbarograph  data  during  eclipses  [see,  e.g.,  Anderson  et 
al.,  1972;  Seykora  et  al.,  1985;  Jones  et  al.,  1992;  Aplin  and 
Harrison,  2003],  Thus  our  result  is  broadly  consistent  with 
observations,  but  seems  to  be  inconsistent  with  standard 
stratospheric  bow  wave  theory. 

[68]  The  observations  of  Anderson  et  al.  [1972]  led 
Chimonas  [1973]  to  the  same  conclusion.  He  proposed  an 
alternative  model  for  these  observations  in  terms  of  V- 
shaped  Lamb  wave  patterns  forced  in  the  troposphere. 
Lamb  waves  are  vertically  external  modes  that  propagate 
laterally  at  the  speed  of  sound.  Lamb  wave  energies  and 
pressure  amplitudes  decay  with  height  above  the  surface, 
and  can  be  efficiently  forced  by  diabatic  forcing  near  the 
surface  [Lindzen  and  Blake,  1972]. 

[69]  NOGAPS-ALPHA  produced  a  strong  surface  air 
temperature  decrease  over  land  (Figure  14),  which  is  also 
broadly  consistent  with  previous  eclipse  observations  and 
relatively  well  understood  in  terms  of  modified  surface  and 
turbulent  heat  transport  [e.g.,  Segal  et  al.,  1996],  This 
temperature  drop  should  also  modify  pressure,  and  indeed 
the  largest  pressure  increases  in  Figure  15  are  also  observed 
over  southern  Africa.  Similar  pressure  increases  over  land 
were  simulated  by  Prenosil  [2000]  using  a  mesoscale  model 
for  the  August  1999  solar  eclipse  over  Europe.  However, 
since  the  temperature  responses  in  Figure  14  are  confined 
to  land  rather  than  moving  with  the  eclipse  shadow,  it  is 
not  clear  whether  this  land-locked  thermal  forcing  could 
generate  the  large-scale  bow-like  pressure  disturbance  in 
Figure  15. 

[70]  Chimonas  [1973]  proposed  Lamb-wave  forcing  due 
to  reduced  shortwave  heating  of  tropospheric  cloud  layers 
within  eclipse  shadows.  Figure  8  reveals  net  reductions 
in  radiative  heating  in  the  troposphere  in  the  NOGAPS- 
ALPHA  simulations  that  are  coherent  over  and  move  with 
the  eclipse  shadow  over  both  land  and  ocean.  Thus  they  are 
a  viable  tropospheric  forcing  term  for  the  large-scale  surface 
pressure  response  in  Figure  15. 

6.  Conclusions 

[71]  Our  NOGAPS-ALPHA  simulations  reveal  clear  evi¬ 
dence  of  a  three-dimensional  bow  wave  response  to  eclipse- 
induced  radiative  cooling  of  the  stratosphere,  as  originally 
predicted  by  Chimonas  [1970]  and  Fritts  and  Luo  [1993]. 
Within  the  eclipse  shadow,  NOGAPS-ALPHA  simulates  an 
induced  radiative  cooling  rate  in  the  stratosphere  that  peaks 
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at  ~27  K  day-1,  a  forcing  2-3  times  larger  than  assumed  in 
these  earlier  stratospheric  bow  wave  models. 

[72]  The  bow  wave  generated  by  this  thermal  stratospheric 
eclipse  forcing  appears  most  clearly  in  NOGAPS-ALPHA 
fields  at  mesospheric  heights  near  the  end  of  the  eclipse  at 
0900-1000  UTC,  a  finding  consistent  with  finite  vertical 
propagation  times  of  gravity  wave  modes  from  their  upper 
stratospheric  source.  These  eclipse-induced  bow-wave 
fields  have  small  amplitudes,  yielding  temperature  changes 
<  I  K  and  horizontal  wind  changes  ^2-3  m  s-1.  These 
findings  support  the  small  middle  atmospheric  responses 
predicted  by  Fritts  and  Luo  [1993]  but  contradict  the  much 
larger  responses  inferred  from  rocketsonde  measurements 
by  Quiroz  and  Henry  [1973]  and  Schmidlin  and  Olsen 
[1984], 

[73]  Atmospheric  surface  pressures  in  NOGAPS-ALPHA 
at  1000  UTC  also  show  a  clear  three-dimensional  bow- wave 
response  of  amplitude  ~0. 1—0.5  hPa,  values  consistent  with 
previous  surface  pressure  measurements  during  eclipse  pas¬ 
sages.  Surface  air  temperatures  in  NOGAPS-ALPHA  show 
drops  of  up  to  4  K  during  the  eclipse  over  land,  and  net 
radiative  cooling  occurs  throughout  the  troposphere  within 
eclipse  shadows.  These  or  other  diabatic  changes  to  the 
troposphere  likely  provide  the  forcing  for  this  wavelike 
surface  pressure  disturbance  [e.g.,  Chimonas,  1973], 

[74]  Future  work  could  seek  to  apply  this  global  modeling 
strategy  to  other  eclipse  events  to  see  if  these  findings  are 
general  or  specific  to  this  particular  eclipse.  A  higher- 
altitude  global  model  that  incorporates  realistic  mesospheric 
and  thermospheric  radiative  heating  and  cooling  would 
yield  more  realistic  simulated  middle  atmospheric  responses 
to  the  eclipse.  Additional  high-resolution  middle  atmospheric 
observations  during  eclipse  passages  would  also  be  useful  in 
investigating  the  origin  of  ongoing  large  discrepancies 
between  rocketsonde  data  and  model  predictions. 


Appendix  A:  Specifying  the  Solar  UV  Eclipse 
Shadow  of  4  December  2002 

[75]  The  standard  USNO  eclipse  metric  is  the  eclipse 
magnitude  (or  phase),  EM,  the  fraction  of  the  Sun’s  diameter 
occulted  by  the  lunar  shadow  at  any  given  geographical 
location.  If  the  apparent  diameters  of  the  lunar  and  solar 
disks  are  DL  and  Ds  =  2 Rs,  respectively,  and  the  separation 
between  their  centers  is  6,  then 
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Ds 


(Al) 


[76]  The  maximum  eclipse  magnitude  (EM) max  occurs 
when  <5  =  0  (concentric  lunar  and  solar  disks).  For  total 
eclipses,  ( EM)max  >  1  since  DL  >  Ds:  for  the  4  December 
2002  total  eclipse,  (EM) max  =  1.02437  [Espenak  and 
Anderson,  2001], 

[77]  Espenak  and  Anderson  [2001,  Figure  1]  plots  paths 
of  constant  EM  in  intervals  of  0.2,  rather  than  just  the 
limiting  umbral  (EM=  1)  and  penumbral  (EM  =  0)  paths 
given  in  Figure  la.  Those  data  show  that,  within  the  shadow 
region,  EM  varies  approximately  linearly  with  great  circle 
distance  from  the  central  total  eclipse  region  to  the  penum¬ 
bral  perimeter,  and  so  we  specify  EM  in  this  way  using 
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Figure  Al.  Geometry  of  the  eclipse  as  a  function  of 
distances  (X,  Y)  from  the  center  of  the  solar  disk,  normalized 
by  solar  disk  radius  Rs,  for  EM  =  0.355  and  (EM) max  = 
1 .02437.  Lunar  disk  (shown  in  gray)  moves  along  X axis  and 
progressively  obscures  the  solar  disk.  Obscured  surface  area 
A0b  is  shaded  dark  gray.  Contours  on  solar  disk  show 
normalized  UV  limb  darkening  contours  from  equation  (All). 
Contour  interval  is  0.05  and  maximum  value  is  0.975. 


linear  interpolation  along  great  circle  paths  from  the  central 
umbral  region  to  the  penumbral  perimeter. 

[78]  This  interpolation  procedure  is  straightforward  when 
all  parts  of  the  eclipse  shadow  strike  the  surface  of  the  Earth 
so  that  the  penumbral  perimeter  is  closed,  as  in  Figure  lg. 
During  the  beginning  and  end  phases  of  the  eclipse, 
however,  only  a  part  of  the  two-dimensional  penumbral 
disc  (defining  where  the  Moon  obscures  the  Sun)  actually 
strikes  the  Earth’s  surface,  which  causes  the  penumbral 
perimeter  to  be  open  as  the  partially  eclipsed  Sun  disappears 
from  view  at  the  solar  tenninator  (see,  e.g.,  Figures  lb-le 
and  1  i— 11).  Furthermore,  at  the  very  early  and  late  stages  of 
the  eclipse,  there  are  no  regions  of  totality  within  the 
penumbral  region,  and  so  the  location  and  magnitude  of 
the  largest  EM  value  are  not  defined. 

[79]  To  estimate  EM  values  within  these  regions,  we  apply 
the  following  algorithm.  First,  we  locate  the  open  end  points 
of  the  penumbral  perimeter  at  the  terminator,  compute  the 
great  circle  distance  and  bearing  angle  between  these  points, 
and  use  those  values  to  locate  the  midpoint.  If  there  is  no 
region  of  totality,  we  use  this  midpoint  as  the  geographical 
location  of  maximum  EM,  and  set  its  value  by  interpolating 
linearly  in  time  between  the  nearest  time  of  total  eclipse 
( EM  =  1)  and  the  nearest  time  of  zero  eclipse  ( EM  =  0). 
Thus,  at  a  point  in  time  halfway  between  the  times  of  first 
partial  eclipse  and  first  total  eclipse,  we  set  the  maximum 
Em  =  0.5.  To  set  all  the  other  EM  values  on  the  sphere,  we 
first  artificially  close  the  penumbral  perimeter  by  locating 
the  point  on  the  penumbral  perimeter  farthest  away  from 
this  midpoint.  Using  the  bearing  angle  between  these  two 
points,  we  locate  another  point  in  the  opposite  bearing 
direction  located  at  a  distance  from  this  midpoint  equal  to 
the  separation  between  the  penumbral  endpoints.  This 
projected  point  always  lies  on  the  sphere  within  regions  of 
negative  cos  <j>  at  an  orientation  roughly  equal  to  the 
elliptical  long  axis  of  the  eclipse  shadow  region.  We  then 
form  lines  from  this  projected  point  to  both  penumbral 
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Figure  A2.  Plots  as  a  function  of  Em  of  eclipse  magnitude 
Em  (dotted  line),  obscuration  Eob  (dashed  line),  and  solar 
UV  obscuration  Euv  (thick  solid  line),  for  (EM) max  = 
1.02437.  See  text  for  details. 


endpoints  using  great  circle  projections,  thereby  creating  an 
artificially  closed  penumbral  perimeter.  We  then  set  all  Em 
values  within  this  (now  closed)  perimeter  using  linear 
interpolation  between  the  maximum  Em  value  (either  at 
the  umbral  location  or  at  the  midpoint)  and  all  points  on  this 
perimeter. 

[so]  With  Em  values  on  the  sphere  now  set,  we  use  them  to 
specify  the  reduction  in  solar  UV  irradiance.  For  this,  the 
more  pertinent  parameter  is  the  fraction  of  the  surface  area  of 
the  solar  disk  obscured  by  the  Moon,  known  as  the  eclipse 
obscuration  [Espenak  and  Anderson,  2001],  and  given  by 


Eob  — 


Aob 

*RY 


(A2) 


where  A  ob  is  the  surface  area  of  the  solar  disk  obscured  by 
the  Moon,  shaded  dark  gray  in  Figure  A I .  For  0  <  EM  <  1, 
trigonometric  manipulation  using  the  geometry  in  Figure  A1 
yields  the  following  analytical  expressions: 

-nEob  =  2(jg)  <P  +  -  X  sin  if’,  (A3) 


X  -  1  +  -  2Em  -  2(EM)m lax-  2 EM, 


(A4) 


cj>  =  2  arcsin 


(A5) 


if)  =  2  arcsin 


(A6) 


The  dotted  and  dashed  lines  in  Figure  A2  compare  EM  and 
Eob,  respectively,  for  (EM) max  =  1.02437. 

[8i]  If  solar  irradiance  was  uniform  across  the  solar  disk, 
Eob  could  be  used  to  scale  down  the  solar  constant  So  in  the 
model  within  the  eclipse  shadow  as 


S  =  S0{l-EOb).  (A7) 


While  this  is  a  fair  approximation  in  the  near-infrared,  at 
shorter  wavelengths  solar  intensities  are  brighter  at  the  disk 
center  and  dimmer  near  the  limb.  Furthermore,  the 
magnitude  of  this  systematic  center-to-limb  decrease  in 
solar  intensities  increases  as  wavelengths  decrease  down 
through  the  visible  and  into  the  mid-UV  [e.g.,  Neckel, 
2005].  Given  this  wavelength  dependence,  to  simplify 
things  we  assume  here  that  the  dominant  radiative  heating 
influence  occurs  in  the  middle  atmosphere  at  solar 
wavelengths  200-320  nm  because  of  ozone  absorption  in 
the  Hartley  and  Huggins  bands. 

[82]  Observational  specifications  of  solar  limb  darkening 
generally  come  as  profiles  of  center-to-limb  intensity  var¬ 
iations  of  the  normalized  form 

(A8) 


where  /(A,  p)  is  the  solar  intensity  at  wavelength  A  at  a 
position  p  =  cos  9  on  the  solar  disk,  and  9  is  the  heliocentric 
angle:  this  coordinate  is  related  to  the  radial  distance  R  from 
the  center  of  the  solar  disk  as  R  =  Rs  sin  9.  Here,  we  average 


Solar  Limb  Darkening:  210-330  nm 


0.0 

0.0 


■  Analytical  Form  Used  Here: 

■  Bonnet  (1968):  210-284  nm 

-  Neckel  and  Labs  (1994):  300-330  nm 

Koepke  et  al.  (2001):  250  nm 

-  Hestroffer  and  Magnan  (1998):  303-320  nm 

Paltridge  and  Platt  (1976):  300  nm 
•  Kjeldseth  Moe  and  Milone  (1978):  210-320  nm 

0.2  0.4  0.6  0.8 

R/Rs  =  sin  0  =  (1  -|x2),/2 


Figure  A3.  Solar  limb  darkening  profiles  /A(//),  averaged 
over  the  wavelength  ranges  A  indicated,  versus  normalized 
solar  disk  radius  (1  —  /i2)1/2,  from  the  following  sources: 
black  dot-dashed  curve,  observations  from  Bonnet  [1968, 
Table  II];  black  long-dashed  curve,  fits  to  observations  from 
Neckel  and  Labs  [1994,  Table  1];  gray  dot-dashed  curve, 
Koepke  et  al.  [2001,  equations  (2.1)  and  (2.2)];  gray  solid 
curve,  fits  to  observations  using  Hestroffer  and  Magnan 
[1998,  equations  (1)  and  (5)];  black  short-dashed  curve, 
Paltridge  and  Platt  [1976,  equation  (3.1)];  black  dotted 
curve,  fits  to  observations  from  Kjeldseth  Moe  and  Milone 
[1978,  Table  3].  Thick  black  curve  shows  the  analytical  p1'3 
curve  that  we  adopt  in  this  study. 
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individual  profiles  over  some  band  of  UV  wavelengths  A  = 
[Ai,  A2]  to  yield  a  mean  limb-darkening  profile 

h(p)  m  Aa  Ai  f(A,  p)d\.  (A9) 

[83]  Figure  A3  plots  I\(p)  from  various  published  sour¬ 
ces,  with  our  A  range  for  each  curve  quoted  in  the  bottom 
left  of  the  plot.  Limb  darkening  observations  in  the  200- 
320  nm  band  are  quite  limited.  The  most  complete  measure¬ 
ments  shown  in  Figure  A3  are  those  of  Bonnet  [1968]  and 
Kjeldseth  Moe  and  Milone  [1978],  which  differ  somewhat. 
Greve  and  Neckel  [1996]  reviewed  these  and  other  observa¬ 
tions  at  200-330  nm.  Using  the  ratio  of  the  disk-integrated 
to  disk-center  intensities 


where  u*  =  u(p/p0)y2,  F*  =  F(pl  p0)y2g/(V2d), 

(B2) 

a  =  dz  —  T,  /3  =  dz  +  T,  and  T  =  pj{2p)  +  g/c2  is  Eckart’ s 
coefficient. 

[86]  To  solve  (Bl),  Frills  and  Luo  [1993]  apply  a  Fourier 
transform  in  y  and  a  sine  transform  in  x.  This  leaves  an 
inhomogeneous  ordinary  differential  equation  in  z,  expressed 
by  their  equation  (19),  which  they  solve  numerically  with 
the  aid  of  a  Green’s  function. 

[87]  Flere  we  use  a  computationally  faster  approach  based 
on  a  method  developed  by  Lighthill  [1960,  1978],  A  three- 
dimensional  Fourier  transform  is  applied  to  (Bl)  to  obtain 


F  F{  A) 
Ic  /( A,  1) 


/' 


pH  a,  p)dp, 


(A10) 


Bu*  =  —ik 


(B3) 


Greve  and  Neckel  [1996]  argued  that  the  Kjeldseth  Moe  and 
Milone  [1978]  data  underestimate  the  limb  darkening  (FUC 
too  large),  whereas  the  Bonnet  [1968]  data  fall  toward  the 
low  end  of  the  F/Ic  range  (see  their  Figure  1).  On  the  basis 
of  their  analysis,  we  choose  an  I\(p)  profile  of  the  form 
p}  \ plotted  in  Figure  A3,  that  lies  midway  between  these 
two  limiting  profiles.  This  choice  yields  Fllc  =  0.607,  a 
value  within  the  range  of  the  data  given  by  Greve  and 
Neckel  [1996,  Figure  1].  This  profile  specifies  the  two- 
dimensional  limb  darkening  over  the  solar  disk  as 


/Sz,(7?,<p)=/^{l-[^]2}1/2), 

/  ->\0-65 

=  (l~iR/Rs}2)  ,  (All) 

where  cp  =  arctan (Y/X)  and  R  =  ( X 2  +  Y 2)1/2.  These  limb 
darkening  contours  are  displayed  on  the  solar  disk  in 
Figure  AI. 

[84]  For  a  given  eclipse  magnitude  EM  and  obscuration 
Fob  (A  oh),  we  compute  a  corresponding  obscuration  of  solar 
UV  intensity  as 


j  jAJSD{R,cp)dRdv 

IoS  CiSD(R,<p)dRd<p’ 


(A12) 


Here  u*,  u*  (and  F*,  F*)  are  Fourier-transform  pairs: 


u  (x,y,z)  = 


OO 

///  u* (k,  l,  m)el{la+ly+m:)dk  dl dm.  (B4) 


[88]  The  term  B  in  (B3)  is  obtained  from  L  by  replacing 
dx,  dv,  dz  with  ik,  il,  im,  respectively,  giving 


(B5) 


The  gravity-wave  dispersion  relation  is  B  =  0,  which  can  be 
written  as 


m 


2 


C k2  +  l2)N 2  ^ 

r2  1 


(B6) 


where  u>  =  —kV  is  the  intrinsic  frequency.  Note  that  the 
appearance  of  T2  instead  of  the  more  usual  (pJ2pf  =  T2  + 
N2/c2  in  (B6)  is  a  by-product  of  the  use  of  the  pseudo- 
incompressible  rather  than  the  fully  compressible  equations 
[cf.  Durran,  1989,  equations  (42)  and  (43)]. 

[89]  Dividing  (B3)  by  B  and  substituting  into  (B4)  yields 
the  spatial  solution 


where  the  integration  in  the  numerator  is  performed  over  the 
obscured  surface  area  Aob.  We  evaluate  equation  (A  12) 
numerically  for  the  range  of  eclipse  magnitudes  EM.  The 
resulting  Euv(EfJ)  curve  is  plotted  in  Figure  A3.  We  use  this 
curve  to  convert  our  EM  values  into  Euv  values  for  use  in 
NOGAPS-ALPHA.  These  Euv  values  are  plotted  as  the 
white  contours  in  Figure  1. 


OO 

/// 


—ik  [  a  — 


N2\  F 

B 


e«(fa+fy+mz)  dkdidm.  (B7) 


Poles  in  the  integrand  occur  where  B  =  0,  so  the  integration 
with  respect  to  m  can  be  performed  using  the  method  of 
residues  to  obtain 


Appendix  B:  Fourier  Solution  of  Fritts  and  Luo 
[1993]  Equations 

[85]  From  the  starting  equations  (3)-(7),  Fritts  and  Luo 
[1993]  obtain  relations  for  each  independent  variable,  e.g., 

Lu*m-(a-N2/g)F*,  (Bl) 


PC  /  N^\  Fpimz 

u*  =  2m  /  /  — ikla - J— — e,<'kx+Iy'> dk  dl,  (B8) 

—  OO 

where  Bm  =  2k2m,  and  m  is  now  treated  as  a  function  of  k,  l 
through  (B6).  Lighthill  [1960,  1978]  makes  further  analytical 
progress  by  using  the  stationary  phase  method,  but  this 
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involves  a  far-field  approximation  with  an  unspecified 
region  of  validity  that  is  probably  too  restrictive  for  our 
purposes.  Hence  we  approximate  (B8)  numerically  using  a 
fast  Fourier  transform. 
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